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Abstract 

We present a systematic study of the statistics of the occupation time and 
related random variables for stochastic processes with independent intervals of 
time. According to the nature of the distribution of time intervals, the prob- 
ability density functions of these random variables have very different scalings 
in time. We analyze successively the cases where this distribution is narrow, 
where it is broad with index < 1, and finally where it is broad with index 
1 < < 2. The methods introduced in this work provide a basis for the in- 
vestigation of the statistics of the occupation time of more complex stochastic 
processes (see joint paper by G. De Smedt, C. Godreche, and J.M. Luck 
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1 Introduction 



The question of determining the distribution of the occupation time of simple stochas- 
tic processes has been well debated by probabilists in the past [|l], [|, |], §]. The 
simplest example is that of the binomial random walk, or, in its continuum version, 
of Brownian motion. Let Xt be the position of the walker at time t, where t is an 
integer, for the random walker, or a continuous variable, for Brownian motion. Then 
define the stochastic process o t = signx*. The occupation times T 4 + and Tf are the 
lengths of time spent by the walker, respectively on the right side, or on the left side 
of the origin, up to time t (using continuous time): 

T±=/VI^. (1.1) 



t I — r> 

It is also convenient to define the more symmetrical quantity 

rt 



S t = tM t = f dt'a t , = T+ - T-. (1.2) 



M t , the mean of the stochastic process at, will be hereafter referred to as the mean 
magnetization, by analogy with physical situations where a t is the spin at a given 
point of space. 

A classical result, due to P. Levy, is that the limiting distribution of the fraction 
of time spent by the walker on one side of the origin up to time t, as t — > oo, is given 
by the arcsine lawQ 

lim/,-i T± (x)= 1 (0<x<l) (1.3) 

7T^X(1 - X) 

or equivalently, for the mean magnetization, 

f M (x) = i (-K x < 1). (1.4) 

7TV 1 — X Z 

The interpretation of this result is that, contrarily to intuition, the random walker 
spends most of its time on one side of the origin. Translated into the language of 
a game of chance, this amounts to saying that one of the players is almost always 
winning, or almost always loosing, a situation which can be summarized by the "per- 
sistence of luck, or of bad luck" . 

On the other hand, a well-known result is that the probability for the random 
walker to remain on one side of the origin up to time t, i.e., for the stochastic process 
o t not to change sign up to time t, decays as £~2, which defines the first-passage 



^ee appendix A for the notation. 
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exponent |. The identity between this exponent and that characterizing the singu- 



larity of the probability density functions ( |1.3| , |1.4j) for x — > 1 is not coincidental. 
Both exponents are actually intimately related, as will be explained later in more 
generality. 

The stochastic process a t defined above is simple in two respects. Firstly, the steps 
of the walker, or the increments of Brownian motion, are independent. Secondly, the 
process is zero-dimensional, in the sense that it does not interact with other processes. 

The question of phase persistence for self-similar growing systems, which appeared 
more recently in statistical physics, is a natural generalization of the probabilistic 
problem presented above. The physical systems in consideration are for instance 
breath figures, the pattern formed by growing and coalescing water droplets on a 
plane || [J]; systems of spins quenched from high temperature to zero temperature, 
or more generally to the low-temperature phase ||; or the diffusion field evolving in 
time from a random initial condition || . The question posed is: what fraction of space 
remained in the same phase up to time t, or equivalently, what is the probability for 
a given point of space to remain in this phase. In this context, the stochastic process 
a t is, for instance, the indicator of whether a given point of space is in the dry phase 
for the breath figure experiment, the spin at a given site of a lattice for the zero- 
temperature kinetic Ising model, or the sign of the field at a given point in space for 
the diffusion equation. 

In these cases, the probability for a given point of space to remain in a given 
phase is equal to the probability po(t) that the process a t did not change sign up to 
time t, or persistence probability. For these situations, and similar ones, po(t) decays 
algebraically with an exponent 6, the persistence exponent, which plays the role in 
the present context of the first-passage exponent | of the random walk. 

Phase persistence is far more complex than the simple case of the persistence of 
luck for the random walker, because the physical systems considered are spatially 
extended. Thus, the stochastic process a t , which is defined at a given point of space, 
interacts with similar processes at other points of space. As a consequence, the 
persistence exponent 6 is usually hard to determine analytically. 

Pursuing the parallel with the random walk, one may wonder, for the physical 
systems where phase persistence is observed, what is the behavior of the distribution 



of the occupation time (|1 . 1| ) , or of the mean magnetization ([L2] 



For the zero-temperature Glauber-Ising chain [|10||, the diffusion equation pi, ITT 



the two-dimensional Ising model [TJ3], or for growing surfaces |[L3||, non-trivial U- 



shaped distributions /m are observed. The singularity exponent of these distributions, 
as M t — > ±1, gives the persistence exponent, thus providing a stationary definition of 
persistence, even at finite temperature [0 . The only analytical results at our disposal 
for these examples is for the diffusion equation, using the so-called independent- 
interval approximation [ 10| . However, understanding the origin of the existence of a 
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limiting distribution is easy. In the long-time regime the two-time autocorrelation is 
a function of the ratio of the two times. Therefore the variance of M t remains finite 
(see equation ( |9.3|) below and ref . |10| ) . 

Another case of interest is the voter model. The investigation of the distribution 
of the occupation time for this model was introduced by Cox and Griffeath ]TJ| 



see |15|] for a summary). In contrast with the former examples, persistence is not 



algebraic for this case |L6], [17[], and the distribution ju is slowly peaking with time 
((M t 2 ) ~ 1/ hit). This is akin to what occurs in the Ising model quenched from high 
temperature to the critical point. The explanation is the same in both cases: the 
two-time autocorrelation is the product of a scaling function of the ratio of the two 
times by a prefactor depending on one of the two times. In the case of the Ising 
model, this prefactor is related to the anomalous dimension of the field at criticality 
(see equation ( |10.71 ) below and ref. JTB|). As a consequence, the variance of M t is 
slowly decreasing in time ((M t 2 ) ~ t~ 2 ^' UZc ). Let us recall that the voter model is 
critical ||19|| . The slow peaking of Jm for both models is the signature of criticality. 



To summarize at this point, there is a strong analogy between the persistence of 
luck for the random walk, and phase persistence for self-similar coarsening systems. 
However, for the latter, both the first-passage (or persistence) exponent 9 and the 
limiting distributions lim^oo f t -i T ±(x) and /m(^) are difficult to determine analyti- 
cally. 

Therefore, instead of considering spatially extended coarsening systems, we adopt 
the strategy of investigating simpler (zero-dimensional) stochastic processes, where 
the persistence exponent is a parameter in the definition of the model, yet where the 
distribution of the occupation time, or of the magnetization, is non-trivial. So doing 
we shall gain a better understanding of the nature of the persistent events encoded 
in these distributions, and shall be better prepared to investigate the statistics of the 
occupation time for the more difficult statistical mechanical models. 

The present work is devoted to the study of a first example of this category of 
problems, where the stochastic process at is generated by a renewal process, defined as 
follows. Events occur at the random epochs of time t\, • • • , from some time origin 
t = 0. These events are considered as the zero crossings of the stochastic process 
o t = ±1. We take the origin of time on a zero crossing. This process is known as a 
point process. When the intervals of time between events, T\ = t\, r 2 = t 2 — ti, . . ., 
are independent and identically distributed random variables with density p(r), the 
process thus formed is a renewal process. Hereafter we shall use indifferently the 
denominations: events, zero crossings or renewals. 

In this model the persistence probability Po{t), that is the probability that no 
event occurred up to time t, is simply given by the tail probability: 

rod 

p (t)=V(r>t) = J drp{r). (1.5) 



4 



The approach we use is systematic and applies to any distribution p. In what 
follows p will be either a narrow distribution with finite moments, in which case the 
decay of po(t), as t — > oo, is faster than any power law, or a broad distribution 
characterized by a power-law fall-off with index 9: 

jf 'dr p(r) « (Jff (0<tf<2), (1.6) 

where To is a microscopic time scale. If 9 < 1 all moments of p are divergent, while 
if 1 < 9 < 2, the first moment (r) is finite but higher moments are divergent. In 
Laplace space, where s is conjugate to r, for a narrow distribution we have 

Cp(r)=p(s) = o l-(r)s+ 1 -(r")s 2 + --- (1.7) 

For a broad distribution, Q1.6Q yields 

^ ~o\ l-(r)s + as e (1<0<2), 1 j 

with a = |r(l-0)|r o fl . 

The case 9 — h accounts for Brownian motion. Indeed, as is well known, the 

3 

distribution of first-passage times by the origin behaves at large times as r~2 ; i.e., is 
in the basin of attraction of a Levy law of index | [20]. This can be simply worked 



out in the case of the random walk, where the discreteness of time allows a natural 
regularization of the process at short times. Hence, since we are interested in universal 
asymptotic properties of the distribution of the mean magnetization, the description 
of Brownian motion by a renewal process with distribution of intervals given by (|1.6| ) 
or (|1.8|), with index 9 = |, is faithful. 

Finally, let us give the content of the present work and stress its originality. 



Though the study of renewal processes is classical [^0], [H], |22|], fewer references 
are devoted to renewal processes with broad distributions of intervals. Elements on 
this question can be found in refs. [^, ^3|. The present work provides a systematic 
account of the theory. In particular we perform the scaling analysis of the distributions 
of the various random variables naturally occurring in a renewal process, such as the 
number of events between and t, the epoch of the last event before t, the backward 
and forward times. We also analyze the aging (i.e., non-stationary) properties of the 
distribution of the number of events occurring between two arbitrary instants of time 
and of the two-time correlation function, for the case of a broad distribution p with 
9< 1. 

The study of the limiting distribution of the occupation time of renewal processes 
is the subject of the work of Lamperti ][J. His main result is the expression of this 
distribution for a broad law of intervals p, with 9 < 1 (see equation ( [7.1 0| ) below). 
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This result is recovered by a simple method in ref. [23], for the case where p is a stable 
Levy distribution of intervals, with 9 < 1 . However neither the work of Lamperti, nor 
ref. 



24 



contain the analysis of the scaling of this quantity for the case 1 < 9 < 2. 
The present work fills this gap. 

Finally, the methods used in the present work provide a firm basis to the investi- 
gation of the distribution of the occupation time for the process with time-dependent 
noise considered in [25] . This model is a second example of the category of problems 



mentioned above. There again, the existence of a renewal process in the model plays a 
crucial role. This will be the subject of a joint paper [26]. Both examples, the renewal 



processes considered in the present work, and that just mentioned, are deformations 
of the binomial random walk, or in continuous time, of Brownian motion. However 
they lead to non-trivial limiting distributions of the occupation time and of related 
quantities, as t — > oo. 



2 Observables of interest 

Let us introduce the quantities, the distributions of which will be computed in the 
following sections. 

First, the number of events which occurred between and t, denoted by N t , is the 
random variable for the largest n for which t n < t. The time of occurrence of the last 
event before t, that is of the iV t -th event, is therefore^ 

tN = n h \-t n . 

The backward recurrence time B t is defined as the length of time measured back- 
wards from t to the last event before t, i.e., 

B t — t — tjv, 

while the forward recurrence time (or excess time) E t is the time interval between t 
and the next event, 

Et = tjv+i — t. 

The occupation times T t + and Tf, i.e., the lengths of time spent by the cr-process, 
respectively in the + and — states, up to time t, were defined in the introduction as 

'Jo 2 ' 
hence t = T 4 + + Tf . They are simply related to the sum St by 

S t = C dt'at* = T t + - Tr = 2T+ - t = t - 27T. 
Jo 

2 We drop the time dependence of the random variable when it is in subscript. 
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Assume that cr i=0 = +1. Then 



Tt = Ti + r 3 + • • • + r N 1 
= r 2 + r 4 H h Tjv— 1 + ^ J 

and 

T t + = n + r 3 + ■ ■ ■ + tjv-i + £ t 

?r = r 2 + ^4 h — \-t n 



if JV t = 2k + 1 (i.e., tr t = -1) (2.1) 



if JV t = 2fc (i.e., <7t = +1). (2.2) 



Assume now that a t =o = — 1- Then, with obvious notations, the following relation 
holds 

T±(a i=0 = -1) = T?(a t=0 = +1), (2.3) 

hence 

S t (a t=0 = -1) = -5^=0 = +1). (2.4) 

Finally we shall also be interested in two-time quantities, namely the number of 
zero crossings which occurred between t and t + 1', given by N(t, t + t') = N t+ t> — N t , 
and the two-time autocorrelation of the process a t , defined as C(t,t + t') = (a t a t+t i). 



3 Number of renewals between and t 

The probability distribution of the number of events Nt between and t reads 

p n (t) =V(N t = n)=V(t n <t< t n+l ) = (I(t n < t < t n+1 )) (n > 0), (3.1) 

where I(t n < t < t n +i) = 1 if the event inside the parenthesis occurs, and if not. 
Note that to = 0. The brackets denotes the average over n, r 2 , . . . The case n = is 
accounted for by equation ( |1.5| ). 

Laplace transforming equation (|3.1| ) with respect to t yields 



CPn(t)=p n (s) = { j tn+1 dte~ st ) = (e"^^— ) (n > 0), (3.2) 

and therefore 

p n (s) = p( s ) n) —^- (n>0). (3.3) 

s 

This distribution is normalized since J2^=oPn(s) — 1/s. 

From ( ^.3|) one can easily obtain the moments of in Laplace space. For instance 

C(N t ) = J2np n (s)= n P{ [. (3.4) 



n=l 
oo 



s(l-p(s)) 



C(N?) = £ B ^(.) = M£±*Mi. (3.5) 

* n=l S(l-p(s)) 
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We now discuss the above results according to the nature of the distribution of 
intervals p(r). 

(i) Narrow distributions of intervals 

Expanding ( |3.4| ) and ( |3.5| ) as series in s, and performing the inverse Laplace trans- 
form term by term, yields 

W ^oo (7) + C " 

where the constants C\ and C2 can be expressed in terms of the moments of p. More 
generally, all the cumulants of N t scale as t, as would be the case for the sum of 
t independent random variables, and therefore Nt obeys the central limit theorem 
!lL£2fl. 



(ii) Broad distributions of intervals with index 9 < 1 
Using equation (Oh, (|3.4|) yields 



TO ,« *£(«-)', (3.6) 

while ( |3.5|) yields (N 2 ) ~ t 29 . More generally one expects that the cumulant of order 
k of Nt scales as t ke . Indeed, setting 

Nt - f^) X " 

where we keep the same notation X t for the scaling variable, we obtain, using ( |1.8| ) 
and 

^J^LlM, (3.7) 

the limiting distribution of X t ,&s t — > 00, 

dz 



fx(x) = / —z°- l e z - xz (0 < x < 00). (3.8) 

For small x, expanding the integrand in the right side and folding the contour around 
the negative real axis yields 

fx ^ xZo r(i - 9) ~ r(i - 29) + " ' 

At large x, applying the steepest-descent method, we find the stretched exponential 
fall-off 

i/(i-er 



.1; 



which demonstrates that all the moments of X are finite. 

For 6 = 1/2, the distribution of X is given by a half Gaussian: 

—x 2 /4 

fx(x) = ^—=^ (0<s<oo). 
V 71 " 

(iii) Broad distributions of intervals with index 1 < 9 < 2 
Using equation ( |1.8|) , we obtain 

The fluctuations of the variable N t around its mean tj (r) are therefore no longer 
characterized by a single scale. These fluctuations are very large as we now show. 
We set, still keeping the same notation X t for the scaling variable, 

N ' = ¥) + %{- r(l - e) t)) 1,>Xt - (3 ' u) 

Then, using equations (|1.8| ) and (|3.7|) , we obtain the limiting distribution of X t , as 
t — > oo, 

fx(x) = J^~ ZX+zB - (3-12) 
For large positive values of x, fx falls off exponentially as 



For large negative values of x, linearizing the integrand of equation ( pj.12 ) with respect 
to z e , and folding the contour around the negative real axis, we find 



fx(x) 



\x\ 



X — > — oo P( 



As a consequence, (X t ) is finite, and all higher order moments diverge. Actually 
(X t ) vanishes, as seen from p.9j ) and ( |3 • 1 1| ) , because the difference of exponents 



2-0-1/0 = -(6- l) 2 /9 is negative. One can also check on (|3TT0| ) that (X 4 2 ) is 
divergent, since the difference of exponents 3 — 9 — 2/9 = (9 — 1)(2 — 9)/9 is now 
positive. 

In the limit 9 — >• 2, f x is a Gaussian: 

e -x 2 /4 

fx{x) v^" 

Related considerations can be found in ref. PPf, vol. 2, p. 373. 
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4 Epoch of last renewal 



We begin our study of the distributions of the random variables tjy, B t , E t , T t , and 
St appearing in section § by the example of tN, the epoch of the last renewal before t. 

This quantity is the sum of a random number N t of random variables n, T2, • • ., 
i.e., it is a function of the joint random variables {ti,T2, • • • , tjv}- Now, the latter 
are not independent, although the intervals r 1; r 2 , . . . are, by definition, independent. 
This is due to the very definition of Nt. Indeed if, say, Nt takes the value n, then 
the sum t n — Yli=x T i is constrained to be less than t. Therefore, in particular, each 
individual interval r, is constrained to be less than t. 

These considerations will be now made more precise, by the computation of the 
distribution of t^, below, and by that of the joint probability density function of 
{n, r 2 , . . . , tat} and Nt, in the next section. 

The joint probability distribution ft N ,N of the random variables and N t reads 

ft N ,N(t; y, n) = ^-V (t N < y,N t = n) = (6(y - t N ) I(t n < t < t n+1 )} , 
ay 

from which one deduces the density ft N of tN 

d 00 
ft N (t;y) = -T-V (tN <y) = (S(y-t N )) = ft N ,N(t;y,n). 

d y n=0 



In Laplace space, where s is conjugate to t and u to y, 

ft*r.N{t;y,n) = f tK .N{s;u,n) = (e~ utn / 



C f tN , N {t;y,n) = f tNtN (s;u,n) = (e utn / " +1 dte st ) 
t,y \ Jt„ I 



p(s + u) n - — ^ (n>0). (4.i; 



The distribution (|3.3|) of Nt is recovered by setting u = in (4.1). Summing over n 
gives the distribution of tN, in Laplace space: 

C (e- utN ) = ft N (s;u) = 1 (4.2) 

t ^ I 1 — p(s + u) s 

which is normalized since f tN (s;u = 0) = 1/s. 

The case where the distribution of intervals p(r) is broad, with 9 < 1 (see equation 
( p..8|) ), is of particular interest since it leads to a limiting distribution for the random 
variable t N /t. 

In the long-time scaling regime, where t and tN are both large and comparable, 
or u, s small and comparable, we get 

ft^^^s^is + u)- 9 . 
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This yields, using the method of appendix B, the limiting distribution for the rescaled 
variable t -1 ^, as t — ► oo, 

lim f t - HN (x) = ""—x 8 " 1 (1 - x)- e = p e ,i-e{x) (0 < x < 1), (4.3) 
with x = y/t, and where 

r(a)r(6) 

is the beta distribution on [0, 1]. As a consequence 

(t N ) « 0*. (4.4) 

C — >00 

In the particular case # = |, we have 

lim f t -H N (x) = — 1 (4.5) 
7iya;(l - a;) 

which is the arcsine law on [0, 1]. This is a well-known property of Brownian mo- 
tion ga. 



5 Interdependence of {ti, . . . , Tjy} 

The purpose of this section is to identify the interdependence of the N t first time 
intervals t 2 , . . ., Tiv}, a point both of conceptual and of practical importance (see 
e.g. @). 

In order to do so, we compute the joint probability distribution of these random 
variables and of N t , denoted by 

f{T~l,T2,—,TjS[},N 

(t;T 1 ,T 2 ,...,T n ,n). 

Generalizing the calculations done above, one has, in the Laplace space of all temporal 
variables, 




resulting in 

f{r 1 ,T 2 ,...,T N },N{s; Ui,U 2 ,..-, U n , n) = Yl P( S + U i) ( n ^ °)> 

S i=l 
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where the empty product is equal to 1 for n — 0. 

The marginal distribution (|3.3| ) of N t can be recovered by setting all the Ui = 
in the above expression. By inversion with respect to the variables one gets 

l-p(s) n 

f{n,r 2 ,...,r N },N(s; n,T 2 ,..., T n , n) = Y[ p{n)e STt (5.1) 

s i=i 

and finally 



f{T L ,n,...,T K }M t ^ r i' r 2 ,...,T n ,n) = []J p(n) V(N t - tn = 0) Q(t - t n ) 



\i=i 



where Q(x) is Heaviside step function. This expression clearly exhibits the interde- 
pendence of the random variables {r 1; r 2 , . . . ,t n }. 

Let us investigate the distribution of any one of the T;, say T\. We denote this 
random variable by r t in order to enhance the fact that its distribution is constrained. 
By integration of ( |5.1| ) on r 2 , r 3 , . . . , r n , we obtain the joint distribution of N t and r t 
in Laplace space: 



M*> T > n ) = * [ S(t) (n = 0), 

from which one deduces, by summation upon n, that 

e~ ST 1 — d(s) 

fn(s; r) = p(r) + 5(r) (5.2) 

s s 

This can be finally inverted, yielding 

f Tt (t; t) = p(r)G(t -t) + 5(t) p {t). (5.3) 

This result can be interpreted as follows. The a priori distribution p(r) is unaffected 
as long as r is less than the observation time t. The complementary event, which has 
probability Po(t), corresponds formally to r = 0. 
As a consequence of ( |5.2|) and Q5.3Q we have 



C(rt) = - rdrrp(r)e- = --^, (5.4) 
t s Jo s as 

and 



roc rt roc 

(r t ) = / drrf Tt (t;r) = / drrp(r) = (r) - drrp(r), (5.5) 
Jo Jo Jt 

where the last expression holds when the a priori average (r) is finite. 

Let us discuss the above results according to the nature of the distribution of 

intervals p(r). 
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(i) Narrow distributions of intervals 

If the distribution p is narrow, then ( |5.5| ) shows that (r t ) converges to (r) very 
rapidly. For instance, if p is exponential, then (r t ) = [1 — (1 + At)e~ A *]/A: the decay 
of (r) — (r t ), with (r) = 1/A, is exponential. 

(ii) Broad distributions of intervals with index 9 < 1 

If the distribution p is broad, with 9 < 1, then from (|5.4| ) 

9t° 

The interpretation of this last result is that tjq ~ t is the sum of iV f ~ t 61 time intervals 
n ~ t 1 - 9 . However (%) = (N t r t ) ± (N t ) (r t ) (see equations and (pUp). 

(iii) Broad distributions of intervals with index 1 < 9 < 2 
We now obtain 

(rt) « {r)-£Lt-V-*\ 

t^oo 9—1 

showing that (r t ) converges to (r) very slowly. 

6 Backward and forward recurrence times 

The distributions of B t =t—tjj and E t = t/v+i —t can be obtained using the methods 
of the previous sections. We have 

fB,N(t;y,n) = (8(y-t + t n )I(t n <t<t n+1 )), 
fE,N(t;y,n) = (S(y - t n+x + 1) I(t n < t < t n +i)) ■ 

In Laplace space, where s is conjugate to t and u to y, 

C f B ,N(t;y } n) = f BjN (s; u, n) = ( dte~ st e^'-A , 
t,y \Jt n I 

C f E ,N(t;y,n) = f E , N (s;u,n) = ( f n+1 dte~ st e- u ^-A , 
t,y \Jt n I 



hence, for n > 0, 



fB,N{s;u,n) = p(s) ■ . 

s + u 

, p(s)-p{u) 
u — s 



f E ,N{s;u,n) = p(s)'' 



and therefore 



f B {s;u) = — -, (6.1) 

s + u 1 — p{s) 

Je{s;u) = ^-r. 6.2 

s — u 1 — p(s) 
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Equations ( |6.1|) and ( |4.2|) are related by f B (s;u) = ft N (s + u; —u), expressing the 
identity t N + B t = t. 

We now discuss the above results according to the nature of the distribution p(r). 

(i) Narrow distributions of intervals 

For distributions with finite moments, equilibrium is attained at long times, for 
both the backward and forward recurrence times, with a common distribution given 
in Laplace space by 

fB,eq(u) = fE,eq(u) = Hmsf B (s] u) = limsf E (s] U) = 1 f ^ ■ (6.3) 

By inversion we obtain 

/Wv) = hAv) = A r dT P^ = t¥ • ^ 6 - 4 ) 

v) J y V) 



In the particular case where p is exponential, inversion of ( |6.1|) yields, for finite t, 

f B (t; y) = Ae-^6(t - y) + e- xt 5(t - y). 

The weight of the second term is simply V(B t = t) = po(t). Similarly, by inversion of 
(fp.21), we have 

fE(t;y) = fE, cq (y) = fs^y) = Xe~ Xy = p(y). 

(ii) Broad distributions of intervals with index 9 < 1 

We obtain, in the long-time scaling regime, i.e., for u, s small and comparable, 

f B (s;u)*s- 9 ( S + u) -\ 

which, by inversion, using the method of appendix B, yields the limiting distribution 

lim f t -i B (x) = S —x-\l - xf- 1 = (0 < x < 1). (6.5) 

This result is consistent with (|4.3[) . Similarly, 



f E (s; u)^ — 



s a (u — Sj 

yielding 



sin7r# 1 

km ft-i E \x) = j— — r (0 < x < oo). (6.6) 



Let us point out that the limiting distribution of t/t N+1 , as t — > oo, is given by (|4.3|) . 
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For 9 = |, we obtain 



lim f t -i B {x) = — , (0 < x < 1), 

irjx(l-x) 



which is the arcsine law on [0, 1]. Similarly 

lim f t -i E {x) = 1 (0 < x < oo). 

t-*oo 7r(l + X)^/X 

(iii) Broad distributions of intervals with index 1 < 9 < 2 

The backward and forward recurrence times still admit the common limiting dis- 



tribution (6.4) , which in the present case has the following asymptotic behavior 



/B«(y) = /W») 1 « e ^^) > ( 6 - 7 ) 

characteristic of a broad distribution of index 9—1 < 1. In other words, at equilibrium 
the average backward and forward recurrence times diverge. However their long-time 
behavior can be computed from (|6.1| ) and (|6.2|) , yielding 



r lu\ - 1- P(s)+sdp(s)/ds _ p(s) - 1 + {t)s 

i { t} ~ 3*(l-p( S )) ' t ^ ~ S^l-p(s)) ' 

from which it follows that 

6 
/d\ ^ T j.2-0 /p\ ^ ^0 +2-6 

Wt^oo (2-6)(t) ' { t} t~oo (Q-i)(2-6)(t) ' 



7 Occupation time and mean magnetization 

The central investigation of the present work concerns the determination of the distri- 
butions of the occupation times T t and of the sum St (or of the mean magnetization 
M t ). 

Let us denote by /£± and fg° the probability density functions of these quantities 
for a fixed value of a = a t=0 , and by f T ± and fs the corresponding probability 
density functions after averaging over <Jo = ±1 with equal weights: 

f = \u + + n- 

The symmetry properties ( |2.3| ) and (|2.4|) imply that 

f T+ (t;y) = f T -(t;y), f s (t;y) = f s (t; -y). (7.1) 
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Following the methods used in the previous sections, we define the joint probability 
distribution of T t + and N t at fixed ctq as 



r T l N (t; y, n) = (t+ < y, N t = n) = (6 (y - T+) I (t n < t < t 



dy 
and 

oo 

y ) = j2 y> n ) = ( 5 {y- T t 



n+l. 



n=0 

For <t = +1, Laplace transforming with respect to t and y, using equations (|2.1|) and 
(PD, we find 



f+ +iN (s;u,2k + l) = p K+i {s + u)p K {s 
f+ +N {s;u,2k) = p k {s + u)p k {s) 



s + u 

hence, summing over k, 

?+ / \ fl — p(s + u) A , .l — p(s)\ 1 

y s + u s J t — p{s)p(s + u) 

Using the property Tf = t — T 4 + , and ([7. 1[) , we get 

f T +(s;u) = - (f++(s;u) + f£+(s + u; -it)) = f T +(s + u; -u). 

Similarly, using the property 5 4 = 2T+ — t, and (|7.1|), we get 

fs(s;u) = f s (s; -u) = f T +(s-u;2u). 

(Here it) is the bilateral Laplace transform with respect to y. See appendix A 
for the notation.) 

The final results read 

f ( ) 2s (1 - p(s + u)p(s)) +u(l + p(s + it)) (1 - p(s)) 
JT±{ ' J 2s(s + u)(l-p(s + u) / $(s)) 

~ a (1 - p(a + - m)) + u (p(s + u) - p(s - u)) 
Is[,) { S i-u 2 )(l-p(s + u)p(s-u)) ' 1 } 

Let us discuss these results according to the nature of the distribution of inter- 
vals p(r). 

(i) Narrow distributions of intervals 

If the distribution p is narrow, implying that the correlations between the sign 
process a t at two instants of time are short-ranged (see section ||), it is intuitively 
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clear that St should scale as the sum of t independent random variables and therefore 
obey the central limit theorem. In the context of the present work, we are especially 
interested in large deviations, i.e., rare (persistent) events where S t deviates from its 
mean. In particular, the probability that St is equal to t is identical to po(t), the 
persistence probability. These points are now made more precise. 

First, expanding the right-hand side of ( |7.2| ) to second order in u, and performing 
the inverse Laplace transform, yields 

More generally, it can be checked that all the cumulants of St scale as t. Then, 

f s (t;u) = (e~ uSt ) = e Ks ^ ~ e** (u) , (7.4) 



t— >oo 

where Ks(t; u) is the generating function of cumulants of St, and $(w) is determined 
below. By inversion of ( |7.4j ), 

du 



fs(t;y)~ [p^e^M. 

t^oo J 2l7T 



In the large-deviation regime, i.e., when t and y are simultaneously large and x = y/t 
finite, the saddle-point method yields 

f M (t;x) ~ e~^\ 

where 

S(x) = — min (ux + $(«)) (— 1 < x < 1) 

is the large-deviation function (or entropy) for M t . The functions S(x) and $(it) are 
mutual Legendre transforms: 



Ti(x) + = —ux, u = — — , x 



dx ' du 



Equation ( [7.4|) implies that fs{s; u) is singular for s = <&(u), and therefore, using 
(n, that 

p{<S>{u)+u)p{®{u)-u) = l, (7.5) 

which determines implicitly $(w). For k — > 0, using p(u) ~ 1 — u (t) + |w 2 (r 2 ), 
equation ( |7.5| ) gives 

v ; «-o 2 (r) 

hence 

£(x) « ^ ^x 2 , 

2 ((r 2 ) - (r) 2 ) 
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yielding a Gaussian distribution for S t , thus recovering the central limit theorem for 
this quantity. 

In the particular case of an exponential distribution p(r), all results become ex- 
plicit. We have 

f ( \ 2X + S 

and 



$(w) = VA^M?- A, S(x)=A 1-vW . 



(ii) Broad distributions of intervals with index 9 < 1 

Using equation (|1.8|) , we have, in the long-time scaling regime where u and s are 
small and comparable, 



fs{s] u) 



{s + u) e ~ 1 + (s-u) 9 - 1 
(s + u) e + (s-u) e 



which, by inversion, using the method of appendix B, yields, as t — ■> oo, the limiting 
distribution for the mean magnetization M t , 

2smn9 (1-x 2 )^ 1 

jM\x) = n. (7.6) 

7T (l + x) 28 + (1-X) 29 + 2C0S7T#(1 -£ 2 f 

For x — > 0, the expansion 

2tan(vrg/2) / cos 2 (7rg/2) - \ 
/m(X) " * cos 2 M/2) X + "j 

shows that x = is a minimum of /m(^) for $ < $c> while it is a maximum for > 9 C , 



with # c = cos(vrfl c /2), yielding 9 C = 0.594611 |4 
For x — » ±1, /m(^) diverges as 

/m(x) ^^2- s (l T xr. (7.7) 

x— >±1 7f 

Comparing the amplitude of this power-law divergence to that of the symmetric beta 
distribution over [—1,1] of same index, 

r(0)0r 

we have 



±i /?(x) v ' r(20)r(l -0)' 

The amplitude ratio -8(0) decreases from B(0) = 2 to 5(1) =0; it is equal to 1 for 
— | (see equation ( |7.9|) ). 



Equation (|B.5|) yields the moments of /, 



M, 

M 2 ) = l-#, (M 4 ) = 1- 9 -(A-6 2 ), (M 6 ) = 1-A(23- 10# 2 + 2# 4 ), (7.8) 
and so on. 

In the particular case — \i equation ([7.6|) simplifies to 

/m(x)= 1 g = /9(x), (7.9) 
7tv 1 — ar 

which is the arcsine law on [—1, 1]. 

The corresponding distribution for the occupation time is 

lim WOr) = ^ ff , (7.10) 

7T x 2e + (l-x) 29 + 2cos7r^a; e (l-a;) 

a result originally found by Lamperti ||. If 6 1 = |, the arcsine law on [0, 1] is recovered: 

lim f t -i T +(x) = — , 

TTy x(l - x) 

A last point is that, in the persistence region, i.e., for M t — > 1 and t — > oo, 
has a scaling form, as we now show. Two limiting behaviors are already known: 

(i) for M t = 1, f M (t; x) = \vo(f)5(x - 1), with p (t) ~ (t/r )- e ; 

(ii) for t = oo, the limiting distribution /m(^) is given by equation (|7.7|) . 

In order to interpolate between these two behaviors, we assume the scaling form 



ft\ l - e ( t 
f M (t;x)~l — ) h[z = (l-x) — 



in the persistence region defined above, with z fixed. The known limiting behaviors 
imply 

The scaling function h(z) can be computed in Laplace space as follows. We have 

f s (t;u) = (e' uSt ) = (T^e-^ J™dz e^ z h(z). 

Laplace transforming with respect to t, yields, in the limit s + u — > 0, with u = 0(1), 
since z is finite, 

fs(s; u) » a(s + u) 9 ' 1 / dz e UT ° z h(z). 
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By identification of this expression with the corresponding estimate obtained in the 
same regime from (|7.2| ), we have, with v = —2u, 



2(1 -p(v)) Jo V ; V 2 J 

As a consequence, the function h(z) is non universal since it depends on the details 
of the function p(r). Universality is restored only in the two limits considered above 
(see (OH))- 

(iii) Broad distributions of intervals with index 1 < < 2 
Using equation (|1.8|), (|7.2[ ) yields 

2(r) -af(s + M) 9 - 1 + ( S - u) e - 1 ) 

^ M ) w 9/ \ ^ I ^T7 ( 7 - 12 ) 

2 (r) s — a ((s + w) y + (s — w) y ) 

in the scaling regime where both Laplace variables s and w are simultaneously small 
and comparable. Expanding this expression as a Taylor series in u, and performing 
the inverse Laplace transform term by term, we obtain, for the first even moments of 
the random variable St, 

2t 9 



t^oo 



(2-0) (3-0) (r) 



^oo (4-0) (5-0) (t/~°- (7,13) 

This demonstrates that, in the long-time regime, the asymptotic distribution of S t is 
broad, with slowly decaying tails. This distribution can be evaluated as follows. 

We first notice that equation ( pM2j ) further simplifies for \u\ ^> s. In order for 
both variables of fs(s;u) to stay in the appropriate domains, we consider s > and 
u —\u, with uj real. We thus obtain, in the relevant regime (s <C |a>| <C 1), 

fs(s;u) 



with 



i i to > 
s + cwr 



COS(7T0/2) 



(r) v ' y 2r(0)sin(7T0/2) (r) ' 
In other words, we have the scaling s ~ \u\ e , or, for the typical value of St, 

(S t ) typ ~ t 1 /*. 

The distribution of S't is then given, for y small and t large, by the double inverse 
Laplace transform 

toft / p^L^ « / ^eMr-clM'*. 

27T1 J-oo 27T S + C 0> J-oo 27T 
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Setting 

S t =(ct)V e X t , (7.14) 

we conclude that the scaling variable X t has a nontrivial even limiting distribution, 
as t — * oo, with argument x = yj (ct) 1 ^ 9 , 

f x (x) = L e (x) = [ + °° p e iax -^ e . (7.15) 

J -oo 27T 

This distribution is the symmetric stable Levy law Lg of index 8. Expanding the 
integrand in the right side as a Taylor series in x yields the convergent series 

_ . , i ~ , r((2fc + i)/e) 2k 

For large values of x, Lg(x) falls off as a power law, 

w r (S + l)sinH/2) 

x— >±oo 7f 

as obtained by expanding e~' a ' as 1 — | ct | in equation (|7.15|) . As a consequence, the 
second moment of this distribution is divergent. 

To summarize, the bulk distribution of St is given by Lg(x), with the scaling 
( [T.14| ), while the moments of St scale as ( |7.13[) . These two kinds of behavior are 
actually related, as shown by the following simple reasoning. Using (|7.14j) , we have 

(S 2 ) ~ t 2 ' 9 (x 2 ) ~ t 2 ' 9 f C dxx 2 Lg{x) ~ t 2 / 9 x 2 c - 9 , (7.16) 

where x c is an estimate for x in the tails. Since x c ~ S c /t 1 ^ 9 with S c ~ t, by definition 
(pi) of ft, we finally obtain (S 2 ) ~ t 3 " 9 as in (|7~lp . Note that, since 2/8 < 3 - 0, 
the typical value of 5* t 2 is much smaller than its average, at large times. A similar 
situation arises in a model considered in ref. [27]. We are indebted to J. P. Bouchaud 
for pointing this to us. 

This argument can be generalized to the calculation of non-integer moments of 
S t , yielding 

(\S t \ p )~t^\ 

with j(p) = p/9 if p < 8, and j(p) = p+1 — 8 if p > 8, a behavior characteristic of a 
bifractal distribution. 

In the limits 8 — > 1 and 8 — > 2, becomes respectively a Cauchy and a Gaussian 
distribution, 

1 , e~ x2 / 4 

Lf)(x) — > — -r, L«(V 



e-i tt(1 +x 2 Y y ' e^2 2^ 
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8 Number of renewals between two arbitrary times 



Consider the number of events N(t, t + t') = N t +t> — N t occurring between t and t + t'. 
The probability distribution of this random variable is denoted by 

p n (t, t + t')=V (N(t, t + t') = n). 

The time of occurrence of the n-th event, counted from time t, is denoted by t' n , 
with, by convention, t' Q = 0. By definition of the forward recurrence time E t , the first 
event after time t occurs at time t[ = E t , when counted from time t. Hence the time 
of occurrence of the last event before t + t', counted from t, reads 

^N(t,t+V) — E t + T 2 -\ h T N (t,t+f)- 

Therefore, 

Pn (t, t + t')=V (t' n <t'< t' n+1 ) = (l(t' n <t'< t' n+1 )) (n > 0). (8.1) 
In particular, for n = 0, we have 

POO 

p {t,t + t') = V(E t >t')= dyf E (t;y), (8.2) 

Jt' 

which is the persistence probability up to time t + t', counted from the waiting time t. 
In Laplace space, where u is conjugate to t! , we thus obtain 

Cp n (t,t + t') 
t> 

Cp (t,t + t') 

where is the Laplace transform of fE(t;y) with respect to y. 

We discuss the results above according to the nature of the distribution of inter- 
vals p(r). 

(i) Narrow distributions of intervals 

As seen in section |], the renewal process reaches equilibrium at long times, with, 
for the forward recurrence time, 

fE,eq(u) = 1 , fE, m (y) — -j—c [ dT P (r) . 

(r) u (r) J y 

Therefore the equilibrium distribution of the random variable N(t, t + t') no longer 
depends on t. In particular, the average of N(t,t + t') is equal to t'/ (r). 



Pn(t,u) = f E (t;u) p(u) r 
l-f E {t;u) 



, x l - P \U) 



(n > 1), 



Po(t,u) 



1.3) 
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For an exponential distribution of time intervals, we have 
Pn (t,t + t')=e- xt '^- (n>0), 



TV. 



which is independent of t, showing that the Poisson point process is at equilibrium 
at all times. 

(ii) Broad distributions of intervals with index 9 < 1 

We restrict the discussion to the probability po{t, t+t'). In the scaling regime where 
t and t' are large and comparable, we have, according to equations ( |3.2|) and fl6.6|) , 

p [t,t + t) « / dx hm / t -i B (x) = / ax-—— — r 

Jt'/t 7T Jt'/t X (l + X) 

t/(t+t') / t \ 

dx(3 e ,i-e{x) = g x — — . (8.4) 



t + t' 

In particular, in the regime of large separations between t and t + t', we obtain the 
aging form of the persistence probability: 

sin ir6 ft'\ 6 
Po {t,t + t') w — — - . 

1«(«C 7T0 \t J 

For example, for 9 = |, 



/, 2 t2.lt 
Po(t, t + t = — arctan W — = — arcsm W . (8.5) 

One can also compute 

p{s) - p(w) 



£ (N(t } t + t')) =Y,np n {s, 



u 



t ,t> ' ' v ' «( M -s)(i-^H)(i-^(s))' 

or alternatively use (pT6|), to find, in the same regime (1 <f~f), 

, . ... sin7r6> (t + t') e - t e . . 

(N(t, t + t')) « * ^ . (8.6) 



A consequence of (|8.4|) and (|8.6|) is that, for 1 < f < t, the probability of finding an 
event between t and t + t' goes to zero. In other words, in order to have a chance to 
observe a renewal, one has to wait a duration t' of order t. The intuitive explanation 
is that, as t is growing, larger and larger intervals of time r may appear. The density 
of events at large times is therefore decreasing. 
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(iii) Broad distributions of intervals with index 1 < 9 < 2 
First, equations (|6.7|) and (|8.2|) imply that, for 1 < f < t, 

Then, in the scaling regime where t and £' are large and comparable, we obtain, using 
equations (|6.2|) and (|3.3|), 



po(s,u) 



a s - 1 - u 9 - 1 



(t) s(s — u) 
which by inversion yields 



Po(t, t + t')= _ T f (V-0- 1 ) - (t + . (8.8) 



(0 - 1) (r) 

For 1 <C £' <C i we recover (|8.7|) , while for 1 <C i <C t' we obtain 

i«(«c (r) 

Equation fl8.8|) can be rewritten in scaling form as 



"* t+t ')=(^W> t " (, " I, ^(^)' ' 8 ' 9 > 

with 

2 (x) = x 9 - 1 ((1 - x)- {6 - 1] - l) . (8.10) 

9 Two-time autocorrelation function 

The two-time autocorrelation function of the cx-process reads 

oo 

C{t,t + t') = (a t a t+t ,) = £(-l)V(*,0- 

n=0 



In Laplace space, using (|8.3| ), one gets 

CC(t,t + f) = C( t , M )= 1 - /g(t; " ) -/ E ( t ;^ l ~M 



t> ' u ' m(1 + p(w)) 

1 / 2 
- l_/ B t ;u _- — - 

u \ 1 + p(u 



lfl_/ B (t;«)_l_V (9.1) 
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(i) Narrow distributions of intervals 
At equilibrium, using (|6.3|), we have 

u\ (t) u [1 + p{u)) J 
Expanding the right side as a Taylor series in u yields the sum rules 

dt'C eq (t>) = {T \~^ ) \ jf dt't' C eq (t>) = ^ - <^> + (9.2) 

Since the second expression can be either positive or negative, depending on the 
detailed form of the distribution p(r), the equilibrium correlation function C eq (t') is 
neither positive nor monotonic in general. 
For an exponential distribution p, one has 

C(t,t + t') = e~ 2Xt ' 

which is independent of t, reflecting once again the fact that the process is at equi- 
librium at all times. 

(ii) Broad distributions of intervals 

In the regime where s~ii<1, the first term in the first line of equation ( |9.1| ) 
dominates upon the second one for any < 9 < 2, so that we have 

C(t,t + t') wpo(t, t + t'), 



with po(t, t + 1') given by (|J) for 9 < 1, and by (fD| ) or (|J) for 1 < 9 < 2. 

Let us finally remark that, from the knowledge of the two-time autocorrelation 
C(t,t + t'), we can recover the asymptotic behavior of (S 2 ), or of {M 2 ), respectively 



given by equations ( |7.3| ), ( |7.8| ) and ( |7.13| ), according to the nature of the distribution 
of intervals p. Indeed, by definition of S t (see equation ( |1.2| )), we have 

(S?) = fdt 2 fdhCtyxM). 
For a narrow distribution p, this yields, as t — > oo, 

S*)^2t j Q dt'C^t'), 



which, using ( |9.2j ) leads to ([7.3[). For a broad distribution p with 9 < 1, we have, in 
the long-time scaling regime 

[S}) w 2t 2 J o dr 2 £ dn C {^j = t 2 £ dx gi {x), (9.3) 

where gi(x) is given in ( |8.4| ). A simple calculation then leads to (M 2 ) as in ( |7.8| ). For 
a broad distribution p with 1 < 9 < 2, using (p.9|), we find 



l)(3-9) (t) A 
yielding, after integration, the first line of ( |7.13|) 
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10 Summary and final remarks 



This article is devoted to the study of the occupation times T t , and mean magneti- 
zation Mt, of renewal processes, that is, processes with independent intervals of time 
between events, interpreted as the zero crossings of the stochastic process at — ±1. 
We also compute the distributions of the random variables naturally associated to a 
renewal process, such that the number of events N t occurring between and t, the 
epoch tjv of the last event before t, the backward and forward times B t and E t . We 
finally investigate the number of events occurring between two arbitrary instants of 
time, and the two-time autocorrelation function. 

The present work is also instructive for the understanding of the role of correlations 
on the behavior of the distributions of sums of random variables. Here the random 
variables are the signs a t = ±1, and 

S t = tM t = f dt'a t > (10.1) 
J o 

is therefore the sum of temporally correlated random variables. 

Three cases are to be considered for the discussion of the results, according to the 
nature of the distribution p(r) of the intervals of time between the renewal events. 

The case where p is narrow corresponds to the domain of application of the central 
limit theorem. Correlations between values of the sign process a t at two different 
instants of time are short-ranged. All the observables of interest mentioned above 
have narrow distributions, as well. In particular, 

(St)typ ~ t 1/2 , 

and the limiting distribution of S t /t l l 2 is Gaussian. The persistence probability po(t), 
i.e., the probability that no event occurred up to time t, decreases, as t — * oo, faster 
than a power law. 

The case where p is broad, with index 9 < 1, corresponds to a maximum violation 
of the central limit theorem. The particular case where 9 = | accounts for Brownian 
motion. Correlations are long-ranged: 

(<^t+f > ~ 01 (j^) , (10.2) 

where the scaling function behaves as (t'/t)~ e in the regime of large separation (1 <C 
t <Ct'). The number N t of events between and t scales as t e . The random variables 
i/v, B t} E t , T t and S t , have limiting distributions, as t — > oo, once rescaled by t. In 
particular, focusing on St, the scaling behavior (|10.2|) implies that 

(^)typ ~ t. (10.3) 



26 



The limiting distribution of S t /t = M t is given by equation ( |7.6|) . It is a U-shaped 
curve for 9 < 9 C « 0.59, while its maximum is at zero, for 6 > 9 C . It is singular at 
M t = ±1, with an exponent equal to 9 — 1. The probability, Po{t), that no event 
occurred up to time t, decreases as t~ , while the probability, po(t,t + £'), that no 
event occurred between t and t + t', behaves as {t 1 /t)~ e in the regime 1 <C t <C if ■ 
(We remind that in the long-time scaling regime, (a t (Jt+t') ~ Po(t,t + t').) 

The case where p is broad, with index 1 < 9 < 2, is intermediate. Correlations 
are again long-ranged: 

(a t a t+t ,)~r^g 2 (-±-), (10.4) 



t + t', 

where the scaling function behaves as (t' /t)~ e in the regime of large separation (1 <C 
t <C t'). The law of N t , centered at tj (r), and rescaled by t 1 ^ 9 , is broad, with index 
9. As t —>■ oo, B t and E t have equilibrium distributions, which are broad, with index 
9 — 1. The scaling behavior ( |10.4j ) implies that 

Sf) ~ f-\ (10.5) 

while a more complete analysis leads to 

(^)typ ~ t 1 ' 8 , (10.6) 

and shows that the limiting distribution of S t /t l l e is the symmetric stable Levy law 
of index 9. The existence of the two scales (|10.5 ) and ( 10.6 ) is also reflected in 



the bifractality of the distribution of St, in the sense that the non-integer moments 
scale as 

(\S t \*>)~t^\ 

with j(p) = p/9, if p < 9, and ^(p) = p + 1 — 9, if p > 9. Otherwise stated, there is 
no "gap scaling" for the distribution of S t (see below). The behavior of po(t) is the 
same as for the previous case 9 < 1. However, po(t,t + 1'), though still decaying as 
t'~ e in the regime 1 t t', is no longer given by a scaling function of the ratio 



t/t + t', as was the case when 9 < 1 (compare equations (|8.4|) and ( ^.9|) ). 

The three types of behavior summarized above, corresponding to a narrow distri- 
bution p (i.e., formally 9 > 2), to a broad distribution with 9 < 1, and to a broad 
distribution with 1 < 9 < 2, are reminiscent of the three typical behaviors observed 
in the nonequilibrium dynamics of phase ordering, respectively at high temperature, 
at low temperature, and at criticality. Consider for instance a two-dimensional lattice 
of Ising spins at — ±1. At high temperature, correlations between spins are short- 
ranged, and the analysis given above applies, i.e., the central limit theorem is obeyed 
by S t . For a quench from a disordered initial state to a temperature below the critical 



temperature, the autocorrelation has a form similar to ( |10.2| ), namely 

t 



t + t' 
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where m eq is the equilibrium magnetization. The scaling function g 1 * 1 ^ behaves as 
(f /t)- x / z for 1 < t < where 2; = 2 is the growth exponent, and A ~ 1.25 is the 
autocorrelation exponent. Correspondingly ( 10. 3|) holds. The limiting distribution of 



St/t = M t is a U-shaped curve, with singularity exponent, as M t — > ±m eq , equal to 
— 1, where ~ 0.22 is the persistence exponent at low temperature [I2|. 



Finally for a quench at the critical temperature, correlations have a form similar 

to (TO , 

(ww) ~ r^flf* (^) , (10.7) 

where (3 = 1/8 and v = 1 are the usual static critical exponents, and z c « 2.17 is 
the dynamic critical exponent. The scaling function ^2 Smg behaves as {i! /t)~ Xc / Zc for 



1 < t < where A c ~ 1.59 is the critical autocorrelation exponent (see e.g. ref. ||18|| ). 
However, in contrast with the situation encountered in the present work, one expects, 
by analogy with statics, that the distribution of S t is entirely described by one scale, 
given by (St)typ = (S"t) 2 ~ f}-fil vz c. T n this sense, this distribution is a monofractal, 
and gap scaling holds: all moments scale as (\S t \ p ) ~ (<St)typ ~ 

t^ p \ with 7(p) = 

p(l — (3/vz c ). These expectations are confirmed by numerical computations |28| . 

Pursuing the comparison between the two situations, persistence for coarsening 
systems decays faster than a power law at high temperature, and as a power law at 
low temperature, as is the case for a renewal process, respectively for p narrow, and 
for p broad, with 9 < 1. Again the status of the intermediate case is different: while 
for the renewal processes under study Po(t) is still decaying as t~ 9 when 1 < 9 < 2, 
for critical coarsening the decay of the probability that no spin flip occurred up to 
time t, po(t), is no longer algebraic. 

To conclude, as we already mentioned, the methods introduced here provide a basis 
for the study of similar questions for the stochastic process considered in ref. [[25[ , and 
to some extent, for the random acceleration problem [29] . 

After completion of this article, we became aware of recent works in probability 
theory on various quantities derived from Brownian motion, whose distributions are 
proven to be beta laws, generalizing thus Levy's arcsine law |I(J . 
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A A word on notations 



Probability densities 

Consider the time-dependent random variable Y t . We are interested in the dis- 
tribution function of this random variable, V(Y t < y), and in its probability density 
function, denoted by f Y , dropping the time dependence of the random variable when 
it is in subscript, 

f Y (t;y) = -^V(Y t <y). 

Time t appears as a parameter in this function. 
Laplace transforms 

Assume that Y t is positive. We denote the Laplace transform of f Y (t;y) with 
respect to y as 

C f Y (t;y) = f Y (t;u) = (e-^\ = 'dyer»y f Y (t;y), 
y * 'Jo 

and the double Laplace transform of f Y (t; y) with respect to t and y as 

Cf Y (t-y)=cU~ uYt ) = f Y {s-u)= / dte~ st dyeT u * f Y (t;y). 

t,y t \ ' JO J 

In this work we encounter random variables Y t (such as St or M t ) with even 
distributions on the real axis, i.e., f Y (t;y) = f Y (t;—y). For these we define the 
(bilateral) Laplace transform as 

£f Y (t;y) = (e-^) = / dye""" f Y (t;y). 
y x 1 J-oo 

Limiting distributions 

In a number of instances considered in this work, Y t scales asymptotically as t. 
Therefore, it is natural to define the scaling variable X t = Y t /t, with density fx(t] x) = 
f t -i Y (t; x = y/t). As t — > oo, this density converges to a limit, denoted by 

fx(x) = lim f t -i Y (x). 

t — ^oo 

These considerations hold similarly for other scaling forms of Y t . 

B Inversion of the scaling form of a double Laplace 
transform 

Consider the probability density function f Y (t;y) of the random variable Y t , and as- 
sume that its double Laplace transform with respect to t and y, defined in appendix A, 
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has the scaling behavior 



Cfy(t;y) = f Y (s;u) = -g(-) (D.i; 

t,y S \S 



in the regime s,u — > 0, with u/s arbitrary. Then the following properties hold, as 
shown below. 

(i) The random variable X t = Y t /t possesses a limiting distribution when t — > oo, i.e., 

f x (x) = \imf t -i Y (t;x = y/t). (B.2) 

t— >oo 

(ii) The scaling function g is related to fx by 

(hi) This can be inverted as 

fx(x) = -— limlm 5 f i-) . (B.4) 

(iv) Finally the moments of /x can be obtained by expanding g(y) as a Taylor series, 
since ( |B.3[ ) implies that 

oo 

<7(£) = (B.5) 

fc=0 

First, a direct consequence of the scaling form ( |B.1| ) is that Y t scales as t, as 
can be seen by Taylor expanding the right side of this equation, which generates the 
moments of Y t in the Laplace space conjugate to t. Therefore (p.2| ) holds. 

Then, ( p.3| ) is a simple consequence of ( p.2| ), since 

f Y (s;u)= dt e~ st (e~ uY ) = / dte~ st (e~ utx " ' 



s + uX 



Now, 

The right side can be rewritten using (p.3j ), yielding ( p.4j) . 
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